Crystal Field and Dzyaloshinsky-Moriya Interaction in orbitally ordered 

La .95Sro.o5Mn0 3 : An ESR Study 
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We present a comprehensive analysis of Dzyaloshinsky-Moriya interaction and crystal-field param- 
eters using the angular dependence of the paramagnetic resonance shift and linewidth in single crys- 
tals of Lao.gsSro.osMnOs within the orthorhombic Jahn- Teller distorted phase. The Dzyaloshinsky- 
Moriya interaction (~ 1 K) results from the tilting of the MnOe octahedra against each other. The 
crystal-field parameters D and E are found to be of comparable magnitude (~ 1 K) with D ~ —E. 
This indicates a strong mixing of the |3z 2 — r 2 ) and \x 2 — y 2 ) states for the real orbital configuration. 

PACS numbers: 76.30.-v, 71.70.Ej, 75.30.Et, 75.30.Vn 



INTRODUCTION 

The importance of orbital degrees of freedom in under- 
standing the complex phase diagrams of the manganites 
is subject of intense research activities (see e.g. refer- 
ences [|| |[ D for an overview). The antiferromagnetic 
insulator LaMn03 (Tn = 140 K) is an orbitally ordered 
system 0], which has been established experimentally 
by resonant X-ray scattering H and neutron diffraction 
0. Moreover, Saitoh et al. recently reported evidence 
for orbital excitations by Raman spectroscopy The 
orbital order in LaMnC>3 is induced by the cooperative 
Jahn- Teller (JT) effect of the Mn 3+ ions (electronic con- 
figuration 3d 4 : t^gCg, spin S = 2), which at tempera- 
tures T < Tjt = 750 K leads to a strong orthorhombic 
distortion of the perovskite structure. In the paramag- 
netic state electron spin resonance (ESR) reveals a single 
exchange-narrowed resonance line with a g value near 2.0 
due to all Mn 3+ ions || and hence directly probes the 
spin of interest. Doping divalent ions like Sr or Ca onto 
the La 3+ place gradually suppresses the JT distortion 
and leads to a ferromagnetic insulating and finally to a 
metallic phase at approximately 15% Sr doping. 

In a previous work, we presented a systematic ESR 
study in single crystals of Lai-^Sr^MnOa with Sr con- 
centrations < x < 0.2 [jio]]. In the JT distorted phase 
the resonance linewidth AH is strongly enhanced com- 
pared to the undistorted phase, reaching maximum val- 
ues of about Ai/ max sa 2.5 kOe. Similar results were 
reported from polycrystalline Lai-^Ca^MnOs jll| and 
oxygen doped ceramic LaMn03+a Moreover, the 

single crystals exhibit a pronounced anisotropy of the 
resonance linewidth in the Jahn- Teller distorted phase, 
which disappears at temperatures T > Tjt- The pure 
LaMnOs sample turned out to be strongly twinned and 
therefore did not allow a detailed analysis of the angular 



dependence. Instead, the La .95Sr .05MnO3 crystal was 
found to be untwinned and can effectively be treated like 
the mother compound (x — 0), as it is still an antifer- 
romagnetic insulator (Tn = 140 K) showing a similar 
magnetic susceptibility. Only the JT transition is shifted 
to lower temperatures, with Tjt (a; = 0.05) = 600 K. 

In a first approach C3), we ascribed the orientation 
dependence of the ESR linewidth to the influence of 
the Dzyaloshinsky-Moriya (DM) interaction, which arises 
from the tilting of the MnOg octahedra along the anti- 
ferromagnetically coupled b axis, only. However, the or- 
thorhombic distortion of the MnC>6 octahedra itself gives 
rise to a crystal-field (CF) induced line broadening of 
comparable order of magnitude and for a complete 
description one has also to take into account the DM 
interaction via Mn-O-Mn bonds within the ferromagnet- 
ically coupled ac plane, which is smaller than along the 
b axis but not negligible. 

In the present paper, we study the angular depen- 
dence of the resonance linewidth in the paramagnetic 
regime at temperatures T > Tn- The contributions of 
CF and DM interaction consist of the superposition of 
the four non equivalent Mn ions in the orthorhombic 
unit cell. The evaluation is based on the structural data 
for LaMnOs determined from neutron scattering exper- 
iments by Huang et al. The application to the ex- 
perimental angular dependence of the ESR linewidth at 
200K and 300K allows to estimate the microscopic CF 
parameters D and E and the DM vectors for all Mn-Mn 
pairs. Throughout this paper we will use the crystallo- 
graphic notation following Huang et at, where the b axis 
(instead of the former c axis in reference |l0|]) denotes the 
direction perpendicular to the ferromagnetically coupled 
ac planes. 
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EXPERIMENT 

The ESR measurements were performed with a Bruker 
ELEXSYS E500 CW-spectrometer at X-band frequency 
(y w 9.35 GHz), equipped with continuous gas- flow 
cryostats for He (Oxford Instruments) and N2 (Bruker), 
which allow to cover a temperature range between 4.2 
K and 680 K. The ESR spectra record the power P ab- 
sorbed by the sample from the transverse magnetic mi- 
crowave field as a function of the static magnetic field H . 
The signal-to-noise ratio of the spectra is improved by de- 
tecting the derivative dP/dH with lock-in technique. A 
small single crystal of Lao.gsSro.osMnOs (volume 2 mm 3 ) 
was fixed on a suprasil quartz rod with low-temperature 
glue (General Electrics), which allowed the rotation of 
the sample around defined crystaUographic axes. We 
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FIG. 1: ESR spectra of La .95Sro.o5Mn0 3 at 200 K (left hand 
side) and 300 K (right hand side) for the magnetic field H 
applied parallel to the crystaUographic axis (a,b,c). Solid lines 
represent the fits using the Dysonian line shape, equation |l| 

measured the angular dependence of the paramagnetic 
resonance with respect to the orientation of the static 
magnetic field H in the three crystaUographic planes (ab, 
ac, and be) at temperatures 200 K and 300 K, deeply in 
the Jahn- Teller distorted orthorhombic phase. Figure |l| 
illustrates typical ESR spectra for the magnetic field ap- 
plied parallel to the three main axes. In all cases, one ob- 
serves a broad, exchange-narrowed resonance line, which 
is well fitted by a Dysonian line shape (Ml given by 



dP d AH + a(H - H rcs ) AH + a{H + H Ies ) 
dH a dH^H - H res ) 2 + AH 2 + (H + H ies ) 2 + AH 2 * 

(1) 

This is an asymmetric Lorentzian line, which includes 
both absorption and dispersion, where a denotes the 
dispersion-to-absorption (D/A) ratio. As the linewidth 
AH is of the same order of magnitude as the resonance 



field H rcs in the present compounds, eq. ([j]) takes into ac- 
count both circular components of the exciting linearly 
polarized microwave field and therefore also includes the 
resonance at reversed magnetic field —H xes . 

Such asymmetric line shapes are usually observed in 
metals, where the skin effect drives electric and magnetic 
microwave components in the sample out of phase and 
therefore leads to an admixture of dispersion into the 
absorption spectra. For samples small compared to the 
skin depth one expects a symmetric absorption spectrum 
(a = 0), whereas for samples large compared to the skin 
depth absorption and dispersion are of equal strength 
yielding an asymmetric resonance line (a = 1). A second 
reason for the asymmetry, which also occurs in insula- 
tors, arises from the influence of non diagonal elements 
of the dynamic susceptibility: In systems with interac- 
tions of low symmetry and sufficiently broad resonance 
lines (H res s» AH) the line shape shows characteristic 
distortions depending on the frequency and orientation 
of the exciting microwave field Jl5| , where eq. (|l|) yields 
a useful approximation of the spectrum. 

Comparing the spectra at 200 K with those at 300 
K one recognizes an increasing asymmetry with increas- 
ing temperature: The average D/A ratio increases from 
about 0.05 at 200 K to 0.35 at 300 K. This can be under- 
stood in terms of the skin effect due to the increase of the 
conductivity with increasing temperature. The resistiv- 
ity p of the sample under investigation is about 500 ficm 
at 200 K and 10 Ucm at 300 K Jl§. We estimated (l(J 
that for resistivities p < 1 Slcm the skin depth becomes 
smaller than the dimensions of the sample, which are 
about 1 mm. As the skin depth is proportional to the 
square root of the resistivity, the influence of the skin ef- 
fect becomes visible even at 300 K. However at 200 K, it 
can be neglected. The remaining asymmetry, which will 
be discussed in connection with the resonance field be- 
low, has to be ascribed to the non diagonal contributions 
of the dynamic susceptibility. 

Figures ^ and || show the full angular dependence of 
linewidth (upper frame), resonance field (lower frame) 
and D/A ratio (middle frame) for 200 K and 300 K, as 
obtained from the fit with eq. |l|. The linewidth exhibits 
a pronounced anisotropy with respect to the crystaUo- 
graphic b axis and a weaker angular dependence within 
the ac plane, which can be empirically described by a 
cos 2 -law. The D/A ratio depends on the orientation of 
the microwave field, which is applied parallel to the ro- 
tation axis, and even becomes negative at 200 K, if the 
microwave field is perpendicular to the ac plane. This 
cannot be understood in terms of the skin effect alone, 
which always produces a positive D/A ratio. Moreover, 
the influence of the skin effect should be negligible at 
200 K. Hence, the asymmetry has to be ascribed partly 
to the influence of the non diagonal elements of the dy- 
namic susceptibility. The angular dependence of the res- 
onance field exhibits additional 90° modulations. They 
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FIG. 2: Angular dependence of linewidth AH (upper frame), 
D/A ratio (middle frame), and resonance field H res (lower 
frame) in Lao.gsSro.osMnOs for the magnetic field applied 
within the three crystallographic planes at 200 K 



FIG. 3: Angular dependence of linewidth AH (upper frame), 
D/A ratio (middle frame), and resonance field H Te s (lower 
frame) in Lao.gsSro.osMnOs for the magnetic field applied 
within the three crystallographic planes at 300 K 



appear to be strongest for the ab plane, where the D/A 
ratio exhibits the most pronounced angular dependence. 
At the same time the D/A ratio is strongly correlated 
with the resonance linewidth, as they both attain their 
maximum values at the same angle. However, for these 
broad signals, the baseline is not so well defined as in the 
case of narrow lines and a shift in the D/A ratio can be 
compensated by resonance field and baseline. For this 
reason, the angular dependence of the D /A ratio can be 
an artificial effect, which is probably related to the large 
linewidth. Hence, we tried to fit the ESR spectra with 
a constant average D/A ratio for each rotation plane. It 
turned out that the spectra at 200 K can be satisfactorily 
fitted as well. This procedure has no visible influence on 
the linewidth, but the resonance field is clearly changed, 
as can be obtained from Fig. ||, which shows the corrected 
data for 200 K. The amplitude of the anisotropy remains 
approximately unchanged, but the additional 90° modu- 
lations more or less disappear. Turning to the resonance 
field data at 300 K, the situation is more complicated, as 
the influence of the skin effect becomes important, too. 
The amplitude of the D/A-angular dependence increases 



and a fit with a constant D/A ratio is not satisfactory 
anymore. Therefore we omit a correction of the reso- 
nance fields at 300 K. 

THEORETICAL BACKGROUND 

A strongly exchange coupled magnetic system like 
LaMnOa can be described by the following Hamiltonian 

(*<i) » 

where the first term describes the the superexchange 
interaction between two next-neighbor Mn spins Si and 
Sj with coupling constant J. The second term describes 
the Zeeman splitting of the spin states with gyromagnetic 
tensor g within an external magnetic field H, where /xb 
denotes the Bohr magneton. The third term Hint in- 
cludes all interactions, which do not conserve the total 
spin and therefore contribute to the broadening of the 
ESR line as there are CF, DM interaction, anisotropic 
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exchange (AE) interaction, dipole-dipole interaction, and 
hyper-fine interaction. Due to estimation of their relative 
strength jllj] it turned out that CF and DM interaction 
yield by far the largest contribution. 



Interactions 

To derive the appropriate expressions for CF and DM 
interaction, we use the structural parameters of the 
LaMn03-IIa sample (space group Pnma) determined by 
Huang et al. fill , because it shows an antiferromagnetic 
ground state and an ordering temperature of about 140 
K, which is consistent with the magnetic properties of our 
sample j[7). Figure 7(a) in reference shows the crys- 
tallographic structure of LaMnC>3 in the strongly Jahn- 
Teller distorted phase. Due to the tilting of the MnOg 
octahedra one can identify four inequivalent positions of 
manganese ions in the orthorhombic unit cell, which is 
illustrated in Fig. 7(b) of the same reference jig]. 



Crystal Field 

In a local coordinate system, where the axes are di- 
rected along the Mn-0 bonds of the MnC-6 octahedra, 
the spin Hamiltonian is usually written in cartesian spin 
components {S x , S V ,S Z ) with parameters D and E as fl§|| 
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Long, middle and short Mn-0 distances in square are 



Rf = {x--Y^ + { V bY + {zcY 

Rt = ix i)V + ( b -f + (zcf 



(6) 
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Ri = (xay + (yby + (- + zYc 
Here capital X, Z and small letters x, y, z denote the 
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(3) 



structure parameters (after Huang et al. |13) for the 
oxygen ions along b and within the ac plane respectively. 



Dzyaloshinsky-Moriya Interaction 

The Hamiltonian, which describes the antisymmetric 
DM interaction fll9fl, can be written as 



where the z axis is directed along the longest Mn-0 
link close to the ac plane and the y axis is parallel to 
Mn-O(l) bond almost along the b direction. Small or- 
thorhombic distortions of the MnOg octahedron will be 
neglected in this section. Because the linewidth is invari- 
ant under shift transformations, we can redefine TL for 
simplification as 



H CF = DS 2 z +E{Sl~S 2 y )-E(S 2 z +Sl + S 2 y ) 

= D'S 2 z +E'S 2 y (4) 

where D' = D — E and E' = —2E. It is useful to recall 
that the y axis is almost parallel to b. The octahedra 
are rotated with respect to their next neighbor within 
the unit cell. Therefore in the crystallographic system 
(a, 6, c) we have different crystal- field components at the 
four inequivalent Mn places which are listed in table Q, 

The CF components for the Mn ion at (0, 0, ^ ) read 



'Hdm — Gij ■ [Si x Sj], (7) 

with the DM vector Gy = dij ■ [noi x noj] being per- 
pendicular to the plane defined by a Mn ion at site («), 
the bridge ligand O, and a Mn ion at site (J), where n 
are unit vectors along Mn 3+ -0 2 ~ bonds |2(], ^l). The 
intrinsic scalar parameter d^ strongly depends on the 
orbital states and the Mn-O-Mn bridge angle. In the 
case of pure LaMn03 both the tilting and the JT dis- 
tortion of the MnOg octahedra account for the origin of 
antisymmetric contributions to the superexchange inter- 
action between the Mn ions. A necessary condition for 
the existence of DM contributions is the lack of a center of 
inversion between the magnetic ions |l9| . With the apical 
oxygen being shifted away from the [010] axis, there is a 
rather strong DM coupling between the ac planes and a 
smaller coupling within the ac planes. Figure ^ depicts 
all next-neighbor couplings that give rise to a DM in- 
teraction. For the Pnma structure, the components Gg- 
(a = x, y, z) of the DM vectors of all Mn pairs within 
the unit cell are listed in table Q. The index j = 1...4 



Resonance Field 



The resonance field of the strongly exchange narrowed 
ESR line is generally determined by the first moment of 
the spectrum. Here wc choose another approach, which 
also allows to take into account the demagnetization ef- 
fect. Hence, we start from the equation of motion p2fl: 



= ig f i B {H - — (2N ZZ - N xx - N yy )}S+ 



dt a ™ 2 

1 

ih 



-[S+,H int ] + i 9f i B M(N xx - N yy )S- (9) 



with the magnetization M and the demagnetizing fac- 
tors N aa . The circular spin operators are defined as usual 
S ± = S x ± iS y . Here we assume an isotropic g value, be- 
cause the influence of the AE interaction can be neglected 
with respect to the crystal field. The DM interaction 
does not contribute to the resonance shift, because the 
expectation value of the commutator ([S + ,H UM ]) van- 
ishes. After linearization the expectation value of the 
commutator ([S + , i/int]) is given by 



FIG. 4: Next-neighbor bonds of the Mn ions, a, /3, and 7 
denote the cartesian components of the DM vector Gij . 



Mat 1 



denotes the four magnetically inequivalent positions of 
the Mn ion in the unit cell. The index i = 1...6 refers to 
the six next neighbors around each Mn site with number 
j, as it is shown in Fig. ||. The absolute values can be 
expressed via two parameters di (inter ac plane) and d 2 
(intra ac plane) as 
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Here M at = { gf i B ) 2 S{S + l)H/[3k(T - 6 C w)] is the 
magnetization per one Mn site, and 0cw = HI K is 
the paramagnetic Curie- Weiss temperature. The crystal- 
field components D a p refer to the coordinate system, in 
which the external field determines the z axis and are 
related to the components D a p in the crystallographic 
system via the usual transformation rules (see Appendix 
2). The sum (j) is running over the four inequivalent Mn 
positions of the unit cell. If we further neglect the de- 
magnetization effect, which according to [ p2| is relatively 
small in the paramagnetic regime, the ESR frequency 
reads 



(hvf = [ gf M B H- 



Mat 



£0«>+£«>-2£Cj>)lll) 



M a t : 
4ff^B 



[£(£i -D$+2iD$)] 



£(D«-fl«-2»4»)] 



where R m measures the Mn-0 distance along the an- 
tiferromagnetically coupled b direction and Ri and R s 
denote the Mn-0 distances within the ferromagnetically 
coupled ac-plane. A microscopic expression for d\ and 
di is discussed in Appendix 1. 



The angular dependence of the resonance field 
H IOS (9,(p) is obtained after substitution of D a by D a p 
(see Appendix 2). Polar angle 9 and azimuth angle ip 
are measured with respect to the c and a axes of the 
orthorhobmic unit cell. 
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Linewidth 



In the case of strong exchange narrowing (Ti. cx 3> 
Hint), the ESR linewidth AH is determined by the sec- 
ond moment Mi of the resonance line divided by the 
exchange frequency u cx p3] 



AH 



1 Mg 



(12) 



In the present case the second moment is determined 
by the sum of DM and CF contributions (M 2 DM + M 2 CF ) 

Within the coordinate system, where the z axis is de- 
termined by the external magnetic field, the second mo- 
ment due to Dzyaloshinsky-Moriya interaction is calcu- 
lated as |E1 M 



M° M = ^5(5 + [(G* ) 2 + (Gij) 2 + 2(G^.) 2 ] (13) 

i,0 



The index j — 1...4 is running over all four magnet- 
ically inequivalent positions of the Mn ion in the unit 
cell. The sum over index i = 1...6 refers to the six next 
Mn neighbors around each Mn site with number j. After 
transformation into the crystallographic system and us- 
ing the values listed in table [n], the average over all four 
positions yields the angular dependence 



M 2 DM = -5(5+ l){(2a 2 +4a 2 h )[l + sin 2 0cos 2 ^] 



+ 4(3 2 ab [1 + sin 2 9 sin 2 <p\ 
+ (2 7fc 2 + 4 7a 2 fc )[l + cos 2 0]} 



(14) 



For the crystal field, the expression for the second mo- 
ment is given in terms of D a p by 



M 2 CF 



45(5+1) -3 
20-4 



3 



Dl 



+ h[(b wx -b m f + wl y ]} 



f)(3)Y2 
VV ' 



(15) 



and again the sum over j is running over the four in- 
equivalent Mn positions in the unit cell. The prefac- 
tors fk allow to separate secular (/1) and non secular 
(/2 and /s) contributions. Again, the angular depen- 
dence is obtained by transforming the CF components 
into the crystallographic system and using the values 
listed in table |. The full expressions divided in secu- 
lar and nonsecular parts are given in Appendix 2. It is 
important to point out here that both secular and non- 
secular parts contain strong contributions proportional 
to cos Aip and sin 40, which should result in a 7r/2 peri- 
odic modulation of the angular dependence. However, in 



the sum (secular) + Mfjr (nonsecular) these contri- 
butions cancel each other for the case /1 = fi = fy. As 
the experimental data do not exhibit any tt/2 periodic 
modulation, only the latter case yields the appropriate 
description of the observed angular dependence. Under 
these circumstances the crystal field contribution reads 
e.g. for the ab plane 

M 2 CF (a6) = ^[45(5 + l)-3]{i^[2^)-^-^] : 



(ih2 



3 

E[3pii ) ) 2 -3(Z?«) 2 



(16) 



(2D<£ - D<£ - D$){dW - D$) cos2^} 



The respective expressions for ac and be plane are ob- 
tained by permutation of {x, y, z) and exchange of tp — ► 8, 
as indicated in Appendix 2. 



ANALYSIS AND DISCUSSION 

As it is shown in Fig. [5| the linewidth data at 300 K 
are well described by the superposition of CF and DM 
contribution to the second moment (Eq. p"2| ). 

The superexchange integral J was estimated in mean- 
field approximation from the Curie- Weiss temperature 
@cw = C(4J ac + 2 Jfc) with ferromagnetic in-plane cou- 
pling J ac to four neighboring Mn ions and antiferromag- 
netic inter-plane coupling J b to two neighboring Mn ions 
and the Curie constant C = 5(5+ l)/3fcB (kn- Boltz- 
mann constant). Assuming — J ac = J b = J and insert- 
ing the Mn 3+ spin 5 = 2 and 6 C w = 111 K @, we 
obtain J = 14 K. Both DM and CF interaction are of 
equal order of magnitude about IK. The DM interaction 
di = 1.00(5) K along the antiferromagnetically coupled 
b axis is about three times larger than d<z = 0.30(3) K 
within the ferromagnetically coupled ac plane. The ab- 
solute values of the crystal field parameters D = 0.73(2) 
K and E = —0.63(2) K are of nearly equal strength indi- 
cating comparable axial (D) and rhombic (E) distortions 
of the MnOg octahedra. 

The angular dependence of the resonance linewidth at 
200 K is well fitted with similar but slightly smaller pa- 
rameters (di = 1.00(5) K, d 1 = 0.26(3) K, D = 0.61(3) 
K, E = -0.58(3) K) in Fig. |. In addition, the reso- 
nance field (g — 1.965) can be described by similar crystal 
field parameters as the linewidth, thereby showing good 
agreement of static and dynamic susceptibility. As the 
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FIG. 5: Angular dependence of the ESR linewidth AH in 
Lao.gsSro.osMnOs for the magnetic field applied within the 
three crystallographic planes at 300 K. The solid lines repre- 
sent the fit with equation The lines below illustrate the 
contributions of CF and DM interaction, respectively 



FIG. 6: Angular dependence of the ESR linewidth AH (up- 
per frame) and corrected resonance field H Yes (lower frame) 
in Lao.gsSro.osMnOs for the magnetic field applied within the 
three crystallographic planes at 200 K. The solid lines repre- 
sent the fit with equations ll| and 111 respectively. 



derivation of the linewidth was carried out in the high- 
temperature limit (T 3> Ocw)j where it approximates 
an asymptotic value Aif(oo), we have to take into ac- 
count, that the fit parameters D, E, d\ and di should re- 
flect the temperature dependence. According to Huber et 
al. [ jnj far apart from any magnetic or structural transi- 
tions, the temperature dependence can be approximated 
by the quotient of the single-ion Curie susceptibility and 
the experimental Curie- Weiss susceptibility as 

AH(T) = T ^ cw Aff(oo). (17) 

Taking into account that the square of CF and DM 
parameters appears in the formulae for the second mo- 
ment, we can extrapolate these parameters to T — ► oo 
via the square root of Eq. |l7|. Using the data at 300 K, 
which is far above magnetic order, we have to multiply 
the parameters by a factor 1.26 and get D(oo) — 0.91 K, 
E(oo) = -0.79 K, di(oo) = 1.26 K and d 2 (oo) = 0.38 K. 
Comparing the values at 200 K, one expects a reduction 
of the parameters by approximately 15% with respect to 



300 K. Indeed, this is well fulfilled, despite the onset of 
the critical behavior on approaching magnetic order at 
140 K, only d\ remains unchanged. 

The equations for the resonance field already contain 
the temperature dependence within the magnetization 
M(T). For this reason it should be better described by 
the high-temperature values of the CF parameters than 
by the temperature dependent ones, which is the case at 
200 K. Due to the rather large uncertainty in the deter- 
mination of the resonance field and the proximity of the 
critical temperature region, this discrepancy is not as- 
sumed to be of significant importance. Thus, at 300 K the 
amplitude of the uncorrected resonance-field anisotropy 
is nicely reproduced by the high-temperature values of 
the CF parameters as £T r ™ ax - H™ n = 95 Oe. 

The strength of CF and DM interaction was in- 
dependently determined from orientation dependent 
magnetization data J26|, measured at 4.2 K in a 
Lao.gsSro.osMnOs single crystal of the same batch as used 
in our ESR experiments. In combination with antiferro- 
magnetic resonance measurements it was shown there, 
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that the ground state (T < Tn) exhibits a canted anti- 
ferromagnetic structure, which rules out any phase sep- 
aration in ferromagnetic and antiferromagnetic regions. 
Within the two-sublattice model Pimenov et al. obtained 
(after translation to our Hamiltonian in units of Kelvin) 
to the following parameters from the magnetization data 
|f§: J = 13.2 K, K z = -2E = 0.93 K, K x = D - E = 
0.95 K, and G = sin(155°)di = 1.18 K. These results 
are in good agreement with our findings, as the values 
of the CF parameters D = 0.48 K and E = -0.47 K 
are essentially equal. The DM contribution is of similar 
strength, where d 2 was neglected. Thus, we could con- 
firm the validity of their model from a microscopic point 
of view. 

ESR studies by Huber et al. in LaMnC>3 Q and Tovar 
et al. ]T3] in the series LaMn03+5 revealed values com- 
parable to our parameters. For pure LaMnC>3, they use 
the CF parameter D — 1.92 K determined by Moussa et 
al. from neutron-scattering experiments |2?j and estimate 
the DM contribution as d w 0.8 K. For an oxygen excess 
6 = 0.03 with a transition temperature Xjt = 600 K 
comparable to Lao.gsSro.osMnOs, the CF parameter D 
is found to be about 20% smaller than in the pure com- 
pound, whereas the DM interaction remains essentially 
unchanged at d 0.8 K. Though these data were mea- 
sured in polycrystalline samples, where no distinction 
between axial(D) and rhombic(i?) CF parameters and 
the inter-plane and intra-plane DM interaction could be 
made, their results agree with ours within a factor of 2. 

The fact, that the two CF parameter are of compara- 
ble strength is of particular importance for the picture 
of orbital order in LaMnC>3. Many illustrations of the 
orbitally ordered state invoke the ordering of d% z 2_ r 2 or- 
bitals (e.g. || ||, ^8|), as originally supposed by Goode- 
nough H in order to describe the antiferromagnetic su- 
perexchange. So far this picture has been widely applied 
and in the absence of an experimentum crucis, which di- 
rectly probes not only the electron-electron interaction 
induced by orbital order (which resonant X-Ray scatter- 
ing is assumed to be able to Q ) but also reveals the kind 
of orbitals involved (which could not be identified by X- 
Ray resonant scattering |2^|), it is the obvious choice in 
view of the observed axial distortion of the MnOe octa- 
hedra. 

Considering the equal contributions of D and E, we 
conclude that this picture of an ordering of d 3z 2_ r 2 or- 
bitals has to be modified as it results in E = in contra- 
diction to our findings, which arc also in agreement with 
the results of Pimenov et al. p3| . In a local coordinate 
system the orbital of the e g electron is a superposition 
4>g.e — Ci0 3z 2_ r 2 ± C2(/) x 2_y2 pCfl , where c\ — 0.8 and 
C2 = 0.6 denote the orbital mixing coefficients down to 
lowest temperatures as reported by Rodriguez-Carvajal 
et al. |jj on the base of neutron-diffraction studies. Mat- 
sumoto p0| found in second order perturbation theory 
the following formula for the CF parameter E 



E=-2VZ— Cl c 2 , (18) 

with the spin-orbit coupling A and the energy separa- 
tion A between the 5 E g ground state and 3 Ti g . This for- 
mula qualitatively confirms our conclusions, as not only 
the sign of our parameters is reproduced but also both 
orbital mixing coefficients must not be negligibly small 
in order to explain the observed rhombic CF parameter 
E. 



CONCLUSION 

We performed a systematic investigation of the angular 
dependence of the paramagnetic resonance in the Jahn- 
Teller distorted orthorhombic phase in Lao.gsSro.osMnOs 
single crystals. We presented a comprehensive analysis 
for the resonance linewidth in high-temperature approx- 
imation, which takes into account the microscopic ge- 
ometry of the four inequivalent Mn positions in the or- 
thorhombic unit cell based on the structural data deter- 
mined for LaMnOs by Huang et al. |l3| . The crystal- field 
parameters for all Mn positions and the Dzyaloshinsky- 
Moriya interaction for nearest-neighbor Mn ions along 
the b axis as well as in the ac plane were success- 
fully extracted as D(oo) = 0.91 K, E(oo) = -0.79 K, 
di(oo) = 1.26 K and d 2 (oo) = 0.38 K. These findings 
shed new light on the microscopic picture of orbital or- 
dering and spin-spin interaction in these compounds. 
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APPENDIX 1: Relative superexchange strength 
and sign of parameters 

We consider the Hamiltonian of superexchange cou- 
pling between two Mn 3+ ions on lattice sites (A) and 
(B) 

H = Jab(Sa ■ S B ) = 2J2 M,Cn( s v ' s c) (19) 
vX 
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where j v ^> ^ v > denote the superexchange parameters via 
one electron states, the symbol r\ refers to the Mn(^4) 
state whereas C refers to Mn(5). In the ground state 
Mn 3+ ions have maximum spin. Therefore we can write: 

J ' 4B= 2^5>«" (20) 

n,i 

Along the b axis, the superexchange is mainly realized 
via the \y 2 — x 2 ) — \y 2 — x 2 ) channel and has antifer- 
romagnetic character. Within the ac plane the ferro- 
magnetic coupling is mainly transferred via the channel 
\y 2 — x 2 ) ~ \y 2 — z 2 ). The relevant states like \y 2 — x 2 ) 
and \y 2 — z 2 ) are mutually orthogonal and yield negative 

Jy 2 —x 2 ,y 2 —z 2 • 

The Dzyaloshinsky-Moriya parameter cIab from equa- 
tion u\ can be derived from the following expression 



J ac ~-^(l){\lcos$ + \ 2 s } 2 (24) 

whereas for DM vector the most important term is 

j y 2_ z 2 y2— z 2 ,y 2 —z 2 ,xy an d therefore 

d2~^^(- A ) 1/2 [Xlcos$ + \ 2 ]\ 7r \ rT (25) 

here # is the angle between uao and Uqb vectors, Jh 
denotes the intra-atomic Hund's exchange parameter, U 
is the charge-transfer energy between the Mn ions. From 
these expression one can clearly see that the relative sign 
of parameters (Lab along the b axis (di in text) and cLab 
within (ac) plane (c?2 in text) are different as well as 
superexchange parameters </& and J ac . 



G 



AB 



2SaSi 



T.hCXr)'(v'\lA~ ) \v) 



E 



r^cc^C'I^IC)} 



(21) 



here i is the imaginary unit, £a and £b are spin-orbit 
coupling constants, A^f denotes the energy splitting be- 
tween the states £ and (,'. I a an d Ib are components of 
the one-electron orbital momentum perpendicular to the 
plane, built up by Mn(A)-bridging oxygen(O) -Mn(B), 
uao and riso are the corresponding unit vectors along 
the Mn-0 bonds. j v ('Xv anc ^ ivCXv' can ^ e expressed via 
a product of dimcnsionless transfer integrals A CT , , A s 
corresponding to Mn-0 bonds |3l]. Q 

In order to explain the main features for the moment, 
we shall simplify the orbital ordering as \y 2 — z 2 ) like 
states along the b axis and alternating ordering like \y 2 — 



2 >, \y 2 



within the (ac) plane. Then one has the 



following situation: 

1) Along the b axis, the largest contribution to the su- 
perexchange comes from j y 2_ z 2y2_ z 2 y 2_ z 2 y 2_ z 2, yield- 
ing 



APPENDIX 2: remarks on the calculation 

Here we outline the transformation between the crys- 
tallographic system [x, y, z) 1 1 (a, b, c) and the coordinate 
system (x,y, z) with the z-axis parallel to the applied 
magnetic field H and rotated by the polar angle 9 and 
the azimuth angle ip with respect to the crystallographic 
system. Assuming an isotropic g value, the DM vector 
transforms like the space coordinates: 



/~<X 



G?- cos 9 cos ip + Gffj cos 9 sin <p — Gy sin ( 
G\, cos ip — G x a sin if 



G'j = G% sin 9 cos ip + G y ij sin 9 sin <p + G z tj cos 6>(26) 

The crystal-field components D a p oc (a(3)/R 2 trans- 
form like the product of the oxygen coordinates a and 
/?. 

For the angular dependence of the resonance field we 
need the following relations: 



1) ab plane 



f) 



J b ~\{\?[\lcosV + \ 2 s ] 2 (22) £(£«+£« - 2 ^) = IJ2(2DW-DW-D$) 

3 ' 

and for the DM vector the most important term is _ ^ _£)(«)) cos 2ip 

jy2- z 2 y 2_ z 2 y 2_ z 2 yz , which leads to 2 ^— ' xx vy 

E(^--^ ±2 ^) = ^E( 2Z ^-^-^) 

rfi~-^( i ) 3/2 ^cos 1 9 + A2]A^A CT (23) 

2) Within the ac plane, the major contribution to the 
superexchange comes from j y 2_ z 2 y 2_ x 2 >y 2_ x 2 y 2_ z 2 and 

hence 2) ac plane (ip = 0) 
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3 i 

i 



i^(I?«- J D«)cos20 



3) be plane = §) 



E(^ + ^- 2 ^) = ^E^-^-^S) 

£p<g-.D«>±2i£<g) = -I^(2£)«-4;)-D«) 

i i 

i 

Concerning the linewidth, the second moment due to 
DM interaction can be expressed in crystallographic co- 
ordinates as 



AL 



DM 



;S(S+l){E(^-) 2 [l + sin 2 0cosV] 



+ ^(Gfjfil + cos 2 0] +E G io G ij sin2^sin 2 
+ E G % G % sin 29cosip + J2 G l G % sin 26> sin ^} 

i,3 i,3 

On evaluation of M^ M , it turns out that all sums over 
cross terms Gf^G^ with a ^ f3 vanish and the remaining 
sums over i, which are equal for all Mn sites j, read: 

E(G^) 2 = 2a 2 + 4a 2 b 

i 

E(^-) 2 =4/3 2 b 

i 

E(G^) 2 = 2 76 2 + 4 7a 2 b (28) 

i 

Concerning the contribution of the crystal field M% , 
it is clear from values listed in Table | that many sums 
over four Mn position vanish in the crystallographic sys- 
tem, which is very convenient. For example 



E 

3 



\-^zz ^yy l-^yz 



= E(^-^)^ 

3 

^ ^ (.Dxx Dyy *)^xy 

3 

= ^2[2D { £ - D ( £ - D$j]D$ 

3 

= ^[2D$ - - Dg]D<£ 

3 



From symmetry point of view it is understandable, be- 
cause there are mirror planes within ab, ac and be planes. 
Having in mind these properties, after some manipulation 
we arrive to the following formulae for the secular part of 
M2 (i-e. when factors fa and fa are equal to zero) and 
nonsecular parts (we consider the case fa — fa): 

1) ab plane (0 = §): 



M 2 CF (sec) 



l[45(5+l)-3]{^[2^-flg-D 



80 
9 



0)1 2 
vv ' 



Et(^-^) 2 + 4 (^) 2 ] 
3 

E[2£« -Dj£- D y ^](D^ - D#)cos2p 
3 

Ytt D ™ - D^f - 4(£«) 2 ] cosM (29) 



Af 2 CF (non) = L [AS {S+l)-3]{\YpD^-D^-D^f 

3 

3 

+ lYttD^f + iD^f] 
3 

+ E^i 2 ^-^-^-]^-^) 



+ 3[p«) 2 -( J D«) 2 ]}cos2^ 
- \YjD^-D^f-±(D<i]f\™^} (30) 



2) ac plane (<p — 0): like ab plane with permutation 
x,y,z — > z, x, ?/ and exchange tp —> 0. 

3) &c plane (ip = -|): like a& plane with permutation 
x, y, z —* z, y, x and exchange ip — > 0. 
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TABLE I: Crystal field parameters for different Mn positions 
in the unit cell. 

D a g\Mn(j)- site (0,0, ±) (±,0,1) (0, ±, ±) (±, ±, 1) 

,(l) 

£,(1) 

D m 

i^xz 

(1) 



D U) 

isxx 


D (l) 


JSxx 


J- 7 XX 


D 


'-'yy 


'-'yy 


'-'yy 


'-'yy 


D 


D (i) 


J- 7 zz 


J-Szz 


iSzz 


D 


D (3) 
J-* xy 


D m 

i-* xy 


D m 

J-* xy 


_£)(!) 
J-* xy 




D U) 

ISXZ 


D m 

-L^XZ 


-D (1) 


D m 

J -^xz 




D U) 

Uyz 


J->yz 


Uyz 


1-1 yz 


D 



TABLE II: DM vector components for different pairs of Mn 
ions in the unit cell of LaMnO;?. 

Mn(i)-site\Mn (j)-site (0, 0, ±) (±, 0, 1) (0, ±, ±) 



j > Vi ~ 2 ' z 
(xj -^,Vi,z. 
(Xj + 2>Vj> z : 



I 

2^ 



(-06,0,-76) (a 6 , 0,-76) (-a 6 ,0, -7i) 

(at, 0,7b) (-a 6 , 0,7ft) (a 6 ,0,7 6 ) 

( — Qac, /3ac, — 7ac) (a ac , — (Sac, 7ac) (ttoc , /3ac , 7ac) 

( — Qac, Pac, 7ac) («oc, ~ /3a c, ~ lac) («ac, /Sac, ™7ac) ( — «oc, — /3ac , 7ac) 

( ftac, /^ac, 7ac) (f^ac, /-^ac, 7ac) (a ac , /^ac, 7ac) ( a ac , ftac> 7ac 

(-Qa C ,/3a C ,7a C ) (Oac, -/Sac, ~7ac) («.c , Pac , ~Jac) (" 



f 1 1 i) 

V 2 ' 2 ' / 

(a 6 , 0, -7 6 ) 
(-a 6 ,0,7 6 ) 



